Simulation of normal distributed smooth fields by Karhunen-Loéve expansion in combination with kriging     

Nils-Otto Kitterr?d  Lars Gottschalk 《Stochastic Hydrology and Hydraulics》1997,11(6):459-482

Simulation of multigaussian stochastic fields can be made after a Karhunen-Loéve expansion of a given covariance function. This method is also called simulation by Empirical Orthogonal Functions. The simulations are made by drawing stochastic coefficients from a random generator. These numbers are multiplied with eigenfunctions and eigenvalues derived from the predefined covariance model. The number of eigenfunctions necessary to reproduce the stochastic process within a predefined variance error, turns out to be a cardinal question. Some ordinary analytical covariance functions are used to evaluate how quickly the series of eigenfunctions can be truncated. This analysis demonstrates extremely quick convergence to 99.5% of total variance for the 2nd order exponential (‘gaussian’) covariance function, while the opposite is true for the 1st order exponential covariance function. Due to these convergence characteristics, the Karhunen-Loéve method is most suitable for simulating smooth fields with ‘gaussian’ shaped covariance functions. Practical applications of Karhunen-Loéve simulations can be improved by spatial interpolation of the eigenfunctions. In this paper, we suggest interpolation by kriging and limits for reproduction of the predefined covariance functions are evaluated.  相似文献   

中国数字地震台网的接收函数及其非线性反演   总被引：36，自引：19，他引：36       下载免费PDF全文

刘启元 Kind  R 《地球物理学报》1997,40(3):356-368

利用中国数字台网（CDSN）记录的85个远震事件的宽频带P波波形数据和分离接收函数的最大或然性反褶方法,获得了CDSN台网10个台站不同方位的岩石层接收函数.利用这些台站不同方位的平均接收函数和非线性接收函数反演方法,获得了上述各台站下方100km深度范围内的岩石层S波速度结构.结果表明CDSN台网各台站下方的地壳厚度和岩石层速度结构存在明显的差别.其中,拉萨台下方的地壳厚度为66kin,壳幔界面较模糊,而余山台下方的地壳厚度为34km,壳慢界面两侧速度反差明显.刮用本文方法估计的地壳厚度与已有的结果基本一致,由于CDSN台网覆盖了中国大陆的各主要构造单元,本文的结果为研究中国大陆的岩石层S波速度结构及其横向变化提供了新的证据.  相似文献   

Determination of S-wave velocity structure using microtremors and spac method applied in Thessaloniki (Greece)     

Paschalis Apostolidis  Dimitrios Raptakis  Zafeiria Roumelioti  Kyriazis Pitilakis 《Soil Dynamics and Earthquake Engineering》2004,24(1):49-67

Array measurements of microtremors at 16 sites in the city of Thessaloniki were performed to estimate the V_s velocity of soil formations for site effect analysis. The spatial autocorrelation method was used to determine phase velocity dispersion curves in the frequency range from 0.8–1.5 to 6–7 Hz. A Rayleigh wave inversion technique (stochastic method) was subsequently applied to determine the Vs profiles at all the examined sites. The determination of V_s profiles reached a depth of 320 m. Comparisons with V_s values from cross-hole tests at the same sites proved the reliability of the SPAC method. The accuracy of the V_s profiles, the ability to reach large penetration depths in densely populated urban areas and its low cost compared to conventional geophysical prospecting, make Mictrotremor Exploration Method very attractive and useful for microzonation and site effects studies. An example of its application for the site characterization in Thessaloniki is presented herein.  相似文献   

Curve-Fitting Schemes for Fragmentation Data     

D. P. Blair 《Fragblast: International Journal for Blasting and Fragmentation》2004,8(3):137-150

Fragmentation measurements in the form of sieve passing and mass fraction data were used to test the capability of three different distributions to fit the observed data over a wide range in fragment size and mass. These distributions were based on Rosin-Rammler, lognormal and simple sigmoidal (S-shaped) functions, having 2 input parameters for the single-component versions and 5 input parameters for the two-component versions. Provided convergence was achieved in the non-linear curve-fitting technique, the two-component versions always provided superior fits to the observed data. However, these versions were very sensitive to variations in the values chosen for the input parameters. In this particular regard, the two-component sigmoidal function was the most robust. The present results also show that the two-component lognormal function provided the best fit to the fragmentation data in a general sense, and the two-component Rosin-Rammler function provided the worst fit. However, there was not a significant difference between any of the three methods.  相似文献   

Seismic studies of the Brasília fold belt at the western border of the São Francisco Craton, Central Brazil, using receiver function, surface-wave dispersion and teleseismic tomography     

Marcelo Assumpo  Meijian An  Marcelo Bianchi  George S.L. Frana  Marcelo Rocha  Jos Roberto Barbosa  Jesús Berrocal 《Tectonophysics》2004,388(1-4):173

The Tocantins Province in Central Brazil is composed of a series of SSW–NNE trending terranes of mainly Proterozoic ages, which stabilized in the Neoproterozoic in the final collision between the Amazon and São Francisco cratons. No previous information on crustal seismic properties was available for this region. Several broadband stations were used to study the regional patterns of crustal and upper mantle structure, extending the results of a recent E–W seismic refraction profile. Receiver functions and surface wave dispersion showed a thin crust (33–37 km) in the Neoproterozoic Magmatic Arc terrane. High average crustal Vp/Vs ratios (1.74–1.76) were consistently observed in this unit. The foreland domain of the Brasília foldbelt, on the other hand, is characterized by thicker crust (42–43 km). Low Vp/Vs ratios (1.70–1.72) were observed in the low-grade foreland fold and thrust zone of the Brasília belt adjacent to the São Francisco craton. Teleseismic P-wave tomography shows that the lithospheric upper mantle has lower velocities beneath the Magmatic Arc and Goiás Massif compared with the foreland zone of the belt and São Francisco craton. The variations in crustal thickness and upper mantle velocities observed with the broadband stations correlate well with the measurements along the seismic refraction profile. The integration of all seismic observations and gravity data indicates a strong lithospheric contrast between the Goiás Massif and the foreland domain of the Brasília belt, whereas little variation was found across the foldbelt/craton surface boundary. These results support the hypothesis that the Brasília foreland domain and the São Francisco craton were part of a larger São Francisco-Congo continental plate in the final collision with the Amazon plate.  相似文献   

航天GPS接收机设计   总被引：1，自引：0，他引：1  

杨春宝  张海云  邱宗德  尤政 《全球定位系统》2004,29(3):19-25

航天GPS接收机为航天器提供航迹、姿态、时问和相对距离等导航信息，提高航天嚣运行的自主性。文中介绍了清华宇航中。航天GPS接收的硬件和软件设计，给出测试结果和分析。  相似文献   

A detailed receiver function image of the sedimentary structure in the Bohai Bay Basin   总被引：5，自引：0，他引：5  

Tianyu Zheng 《Physics of the Earth and Planetary Interiors》2005,152(3):129-143

If the site at which receiver functions are constructed is filled by sediments, then the waveforms from these receiver functions are dominantly controlled by the sedimentary structures within the first few seconds after the direct P arrival. Based on this observation, waveform data collected at 44 temporary seismic stations have been used to image the sedimentary structure of the Bohai Bay Basin, a major continental petroliferous basin in Eastern China. An adapted hybrid global waveform inversion method was applied to the receiver functions to extract structural information beneath each of the stations. The derived S-velocity structure provides for the first time, a basin-scale seismic image of detailed sedimentary stratification. The sedimentary cover of the basin is about 2-12 km thick, consisting of Cenozoic, Mesozoic, and Paleozoic strata from top to bottom. The structural features presented in the S-velocity image coincide quite well with the depression-uplift type of tectonic system in the Bohai Bay Basin. The reconstructed morphology of the sedimentary layers provides seismological evidence for the two-stage evolution of the intracontinental basin that were caused by an intensive tectonic regime transition in late Mesozoic immediately following the lithospheric reforming of the Eastern China continent.  相似文献   

Comment on 'On sensitivity kernels for 'wave-equation' transmission tomography' by de Hoop and van der Hilst     

F. A. Dahlen  Guust Nolet 《Geophysical Journal International》2005,163(3):949-951

  相似文献   

Martian Topography: Scaling, Craters, and High-Order Statistics   总被引：1，自引：0，他引：1  

Vladimir Nikora  Derek Goring 《Mathematical Geology》2005,37(4):337-355

The high-order structure functions of Mars topography reveal three specific ranges of scales: (1) scaling range at small scales where the structure functions exhibit scaling behavior; (2) transition range where the structure functions continue to grow but do not reveal scaling; and (3) saturation range at large scales where the structure functions saturate. The scaling and saturation ranges are explored in detail in respect to scaling and intermittency. Analysis of the Mars Orbiter Laser Altimeter (MOLA) data and computer simulations suggest that there are two potential contributors to the small-scale scaling: (i) scale-invariant surface formation; and (ii) effects of discrete morphological forms such as craters. The crater effect also provides an explanation for the large-scale intermittency revealed using the normalized structure functions within the saturation range, which cannot be explained by the ‘scale-invariant’ concept. Overall, the obtained results suggest that the “crater” contribution to the structure function behavior often dominates over the effect of the scale-invariant surface formation.  相似文献   

Elastic lateral dynamic impedance functions for a rigid cylindrical shell type foundation          下载免费PDF全文

R. He  R.Y.S. Pak  L.Z. Wang 《国际地质力学数值与分析法杂志》2017,41(4):508-526

Elastic lateral dynamic impedance functions are defined as the ratio of the lateral dynamic force/moment to the corresponding lateral displacement/rotation at the top ending of a foundation at very small strains. Elastic lateral dynamic impedance functions have a defining influence on the natural frequencies of offshore wind turbines supported on cylindrical shell type foundations, such as suction caissons, bucket foundations, and monopiles. This paper considers the coupled horizontal and rocking vibration of a cylindrical shell type foundation embedded in a fully saturated poroelastic seabed in contact with a seawater half‐space. The formulation of the coupled seawater–shell–seabed vibration problem is simplified by treating the shell as a rigid one. The rigid shell vibration problem is approached by the integral equation method using ring‐load Green's functions for a layered seawater‐seabed half‐space. By considering the boundary conditions at the shell–soil interface, the shell vibration problem is reduced to Fredholm integral equations. Through an analysis of the corresponding Cauchy singular equations, the intrinsic singular characteristics of the problem are rendered explicit. With the singularities incorporated into the solution representation, an effective numerical method involving Gauss–Chebyshev method is developed for the governing Fredholm equations. Selected numerical results for the dynamic contact load distributions, displacements of the shell, and lateral dynamic impedance functions are examined for different shell length–radius ratio, poroelastic materials, and frequencies of excitation. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献